The Physics of Chromo spheric Plasmas 
ASP Conference Series, Vol. 368, 2007 
Petr Heinzel, Ivan Dorotovic and Robert J. Rutten, eds. 



Observations and Modeling of Line Asymmetries in 
Chromospheric Flares 



Arkadiusz Berlicki 

Astronomical Institute, Wroclaw University, Wroclaw, Poland 

Abstract. For many years various asymmetrical profiles of different spectral 
lines emitted from solar flares have been frequently observed. These asymme- 
tries or line shifts are caused predominantly by vertical mass motions in flaring 
layers and they provide a good diagnostics for plasma flows during solar flares. 
There are many controversial results of observations and theoretical analysis of 
plasma flows in solar chromospheric flares. The main difficulty is the interpre- 
tation of line shifts or asymmetries. For many years, methods based on bisector 
techniques were used but they give a reliable results only for some specific con- 
ditions and in most cases cannot be applied. The most promising approach 
is to use the non-LTE techniques applied for flaring atmosphere. The calcula- 
tion of synthetic line profiles is performed with the radiative transfer techniques 
and the assumed physical conditions correspond to flaring atmosphere. I will 
present an overview of different observations and interpretations of line asym- 
metries in chromospheric flares. I will explain what we have learnt about the 
chromospheric evaporation in the frame of hydrodynamical models as well as 
reconnection models. A critical review will be done on the classical methods 
used to derive Doppler-shifts for optically thick chomospheric lines. In particu- 
lar, details on the new approach for interpreting chromospheric line asymmetries 
based on the non-LTE techniques will be presented. 



1. Introduction 

Spectroscopic observations of solar chromospheric flares show that the line pro- 
files emitted by the flaring plasma almost always exhibit asymmetries or shifts. 
These features are surely due to the chromospheric plasma motion and the re- 
sulting Doppler-shifts effects. Interpretation of the shape of line profiles allows 
us to understand the nature of plasma flows during solar flares. Plasma flows in 
the chromosphere are important in the analysis of dynamics and energetics of so- 
lar flares as well as these phenomena supply matter to the coronal parts of flares 
in the process of chromospheric evaporation. Understanding of the mechanisms 
which generate flows is necessary for complete description of solar flares. 

Spectral line asymmetries arise only because of the line-of-sight component 
of plasma velocity which due to the Doppler-shift is responsible for the modifi- 
cation of spectral line profiles. Therefore, for the flares located on the solar disk 
center we are able to analyze flows oriented perpendicular to the solar surface. 
For simplicity, the theoretical analysis of line asymmetries often assume that the 
emitting region is located on the solar disk center and only vertical flows are 
present. 
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Figure 1. Two examples of the Ha line profiles observed during solar flares. 
Left: shifted profile - the Doppler velocity may be determined from the shift 
of the whole line. Right: complicated asymmetric profile - it is impossible to 
determine the Doppler velocity from such line profile using bisector methods. 



Unfortunately, the interpretation of the shape of line profiles is not trivial 
and to determine the velocity field we must use advanced methods based on 
hydrodynamical and radiative transfer techniques. There are some rare cases 
when the interpretation of line asymmetries is relatively simple. When the 
flow of the whole flaring region is homogeneous then the shape of the line is 
not disturbed compared to the static symmetric line profile, but the whole line 
is just shifted towards longer or shorter wavelengths (Fig. Q] - left). We can 
measure this shift and using a simple Doppler formula calculate specific value of 
the line-of-sight velocity. Unfortunately, in most cases the line profiles of solar 
flares exhibit much more complicated structure (Fig. [1] - right) which suggests 
that the velocity field is not homogeneous and different parts of the flare move 
in different way. 

There are two main approaches to the modeling of chromospheric line asym- 
metries. One is based on the hydrodynamics calculations, where the time evolu- 
tion of the solar flare atmosphere is calculated and the radiative transfer formu- 
lae are used to calculate the radiation from this evolving atmosphere. Another 
approach is based on semiempirical models of solar flares. 

In this paper I will present a short review on observation and interpreta- 
tion of the chromospheric line asymmetries observed during solar flares. The 
term chromospheric lines is commonly used to describe all spectral lines formed 
in the solar chromosphere, where the temperature is around 10 4 K. These lines 
are formed in strong non-LTE conditions and complicated radiative transfer 
calculations are necessary to describe the formation of these lines. Strong chro- 
mospheric lines are usually optically thick what means that the optical thickness 
of the plasma in these lines is very large (r >> 1). 

There is a wide literature concerning this topic and I provide the readers 
some references contained the most important results. I will concentrate on 
chromospheric parts of solar flares where the emission comes from cool (10 4 K) 
plasma. This emission is produced mostly in strong resonance lines of hydro- 
gen, calcium or magnesium (Ha, H/3, H<5, CallH, K, etc.) Since the most 
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spectroscopic observations performed during past years concern Ha line, the 
interpretation of these results will take considerable part of this review. 



2. Early Observations and Interpretations 

Asymmetries of chromospheric lines emitted by the flaring plasma have been 
observed for more than half of the century. After the solar spectrographs were 
developed in order to produce the solar spectrum with sufficient spectral resolu- 
tion, observers noticed that the line profiles emitted during solar flares are not 
symmetric (Fig. [2|). It was clear from the beginning that these kinds of mod- 
ifications of the spectral lines are due to the mass motion driven during solar 
flares. However, the mechanisms which could drive the plasma flows were not 
known yet at those days. 

The first analysis of the line as ymmetries were con centrated on statistical 
description of the behaviour of lines. Svestka et al. ( 19621 ) presented a qualitative 



analysis of 244 Ha and CallK spectra of 92 flares. They found that the blue 
asymmetry (blue wing enhancement) occurs mainly in the early phase of flares, 
before flare maximum. However, only 23% of flares contain at least one region 
with blue asymmetry. 80% of flares exhibit the red asymmetry which dominates 
during and after the maximum of flare. It is worth to notice that only 5% of flares 
shows blue asymmetry exclusively. However, because not all flares were observed 
from their beginning, the occurrence of blue asymmetry may be missed for many 
flares. In Fig. the tim e evolution of the asymmetry is presented. Similar 
analysis was performed by iTangj ( 19831 ). By inspection of off-band filtergrams of 
60 flares obtained in ±1 and ±2 A from the Ha line center he found that 92% 
of flares show red asymmetry and only 5% show blue asymmetry. 

Early statistical analysis gave us an qualitative idea about the line asymme- 
try but the physical interpretation of the shape of line profiles is not considered. 
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Figure 3. Absolute mean value of the asymmetry of the Ha line profiles as a 
function of time development of flares (solid line) . Dashed curve takes positive 
and n egative signs into account. Time of flare maximum = (|Svestka et al.l 

GUI. 



The asymmetric line profile contains an information about the velocity field in 
the region where the observed line is formed. Therefore, the first problem which 
needs to be solved is to determine this velocity field using the observed line pro- 
file. Secondly, we have to answer the question why do we observe such a flow of 
plasma, what can generate the flows? 

Determination of the velocity from the observed profiles of the chromo- 
spheric lines is not a trivial task due to the complicated processes of line forma- 
tion and complex velocity field in the chromosphere. Strong chromospheric lines 
are optically thick and the radiation observed at differe nt parts of the line profile 
comes from different height z across the chromosphere ( Mihalas1ll978l ). The core 
of strong chromospheric lines (hydrogen Balmer lines, Ca lines) is formed much 
higher than the wings of these lines. Therefore, the asymmetry of the specific 
line depends on the relation between the height z of formation of given part of 
the line profile and the value of the velocity at this height. This means that 
if the function describing the velocity across the height in the chromosphere is 
complicated, the emergent line profile has also very complicated shape. 

Despite of all difficulties with determining correctly the velocity from the 
line profile shape, many authors tried to use Ha, Call and other lines to find 
the velocity in the flaring chromosphere. All these determinations were based 
on the measurements of Doppler-shifts of the line cores or, more commonly, on 
the bisector technique. 

Using the shift of the line core to obtain the velocity may be misleading 
and the determined velocit y is not corr ect when the velocity gradient in the 
chromosphere is significant ( Athav 19701 ). Unfortunately, the estimation of the 
Doppler-shift obtained with the bisector technique can also give wrong results. 
The Doppler-shift of the line profile correspond to the shift of the central point 
of the bisector connecting the two wings of the spectral line. Since the different 
parts of the spectral line are formed at different height in the chromosphere, 




Figure 4. Temporal variation of the Ha line profiles observed during solar 
flare. The black arrow indicates the peak of asymme tric lines. The bottom 
profil e represents the emission of the quiet-Sun area (|Ichimoto fc Kurokawal 
Il98l . 



using bisectors connecting the wings observed at different frequency, we can 
estimate the line-of-sight velocity at different layers of the chromosphere. How- 
ever, because the radiation of specific frequency within the spectral line does 
not come from one narrow layer of the chromosphere but rather from geomet- 
rically thick region, we cannot say that the Doppler velocity determined form 
given bisector correspond to plasma flow at given height in the chromosphere. 
Moreover, if the velocity gradient in the chromosphere is large then the bisector 
method cannot be used because the Doppler-shift of any bisector results from a 
superposition of many shifts due to the motion of plasma with different velocity 
along line-of-sight. In spite of this the bisector technique was commonly used 
for many years until it was replaced by more advanced, complicated but much 
more precise non-LTE radiative transfer techniques with velocity field included. 

One of the first interpre tation of chromosp heric line asymmetry observed 
during flare was presented by I Acton et al.l ( 19821 ). The authors postulated that 
these asymmetries are due to ch romospheric eva poration driven by accelerated 
electrons or thermal conduction. Neupert ( 19681 ) was probably the first who re- 
alise that chromospheric plasma heated during solar flare may evaporate. This 
evaporated plasma provid es material for loop prominences often observed as th e 
so-called post-flare loops ( Kopp Sz Pneumanll 19761 ; lAntiochos &: Sturrock| [i978). 
For m ore complete r eview of chromospheric evaporation see Hudson paper in this 
book. Acton et al.l ( 19821 ) postulated that in the analysed flare the non-thermal 
electrons heat the chromosphere mainly during the impulsive phase, while ther- 
mal conduction from the hot coronal plasma heated earlier dominates during 
the late, thermal phase. Both mechanisms drive upflow of the cool plasma. The 
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Figure 5. An example of temporal variation of the downflow velocity ob- 
served in flaring region. Filled circles correspond to the values obtained from 
the shift of the Ha far wings, while crosses to the values obtained from the shift 
of the line peak. Open circles represent the time evolut ion of the Ha intensity. 
The microwave emission at 3750 MHz is also shown (|Ichimoto fc Kurokawal 

QUI). 



authors stress that for the first time they observed chromospheric evaporation 

in Ha line. 

Contradictory results were published bv llchimoto Sz Kurokawa (1984) who 
suggested that during the impulsive phase of solar flares the downflow of the 
cool chromospheric plasma is present. These results are based on the large 
red asymmetry of the Ha line observed during the impulsive phase of many 
solar flares (Fig. S]). The downward motion increases at the onset of a flare to 
its maximum velocity of 40 to 100 km s _1 shortly before the impulsive peak of 
microwave emission, and rapidly decreases before the Ha reaches its maximum 
(Fig. [5]). The red asymmetry of the Ha line may be also explained by the 
attenuation of the blue wing by the rising plasma over the flare but the authors 
exclude this case because the optical thickness of a rising cloud is too small to 
explain the emission deficit of the blue wing of Ha line. Also th e high tempora l 
resolution spectroscopic observations of Ha line performed by IWuelser ( 19871 ) 
confirm the existence of red asymmetry during the impulsive phase of solar 
flares. The largest asymmetry is observed during the maximum of microwave 
emission what confirms the earlier results that the plasma downflow is driven 
by the accelerated non-thermal electrons. 



3. Hydro dynamic Modeling of the Flows 

The downflow of cool plasma in the form of chromosp heric conden s ations ob- 
served during solar flares was predicted theoretically bv lFisher et al. ( 1985al ). If 
a region of the chromosphere heated by non-thermal electrons is thick enough, 
then the rapid temperature increase produces an enhanced pressure in the heated 
region. This overpressure, besides the evaporation, also drives downward moving 
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Figure 6. Time evolution of the temperature and velocity in the loop atmo- 
sphere heated by the low (upper two panels) and the high (lower two panels) 
flux of non-thermal electrons. In the first case the upflow with the low velocity 
in the transition region and in the chromosphere is present (gentle evapora- 
tion). For the strong flux of non-thermal electrons the high- velocity upflow (up 
to — 500 km s -1 ) in the high temperature region is obtained (explosive evap- 
oration) and downflow (up to +40kms~ 1 ) is observed in the c hromosphere 

chrom ospheric condensations). For more details see the paper of lFisher et al.l 

1985ah . 



cool and dense chromospheric condensations ( Fisher et al.lll985bl ) which seem to 
be res ponsible for r e d asym metry of the Ha line profiles reported by many au- 
thors. iFisher et al.1 ( 1985bl ) modeled the hydrodynamic and radiative response 
of the atmosphere to short impulsive injections of non-thermal electron beams 
(Fig. [6]). They showed that a high-energy flux of non-thermal electrons drives 
explosive evaporation accompanied by the formation of cool chromospheric con- 
densations in the flare chromosphere. A different situation occurs when the 
flux associated with non-thermal electrons is very low. Then only a weak chro- 
mospheric evaporat ion takes place. This kind of evaporation is referre d to as 
gentle evaporation ( Antiochos Sturrockl fl978l ; ISchmieder et al.l Il987l ) and it 
can be observed in chromo spheric spectral lines like Ha or in Call 8542 A. 
Antiochos <fe Sturrockl ( 1978! ) suggested that the gentle chromospheric evapora- 
tion may also occur after the primary energy release, when the non-thermal 
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Figure 7. Ha line profiles observed during the gradual phase of a solar 
flare (solid lines). Weak blue-shift of these lines suggests slight upflow of 
the plasma with the velocity of a few kms -1 interpreted as gentle evapora- 
tion. Dashed lines rep resent the reference line profiles of the quiet Sun area 
(|Schmieder et al.lll987l ). 



electron flux is stopped. This evaporation could be driven by the large con- 
ductive heat flux from a high temperature flare plasma contained in magnetic 
tubes above the chromosphere. Such physical conditions may appear during 
the gradual phas e of solar flar e s, wh en there is no significant flux of non-thermal 
electrons. In the lForbes et al.l ( 19891 ) model for flare-loop formation by magnetic 
reconnection the conduction of the thermal energy generated at the slow-mode 
sho cks drives a gentle evap orative upflow from the ribbons. 



Schmieder et al.l ( 19871 ) observed small but long-lasting blue-shifts in flare 



ribbons in the Ha line during the gradual phase of three solar flares and inter- 
preted them as due to upflows with velocities less than 10 kms -1 (Fig. [7]). These 
upflows were believed to be caused by gentle chromospheric evaporation driven 
by the heat conduction along the field lines connecting the chromosphere with 
a reconnection site in the corona. 

The downflow of cool chromospheric plasma during the impulsive phase of 
solar flares predict e d in t he theoretical calculations was reported by many au- 
thors. IZarro et al. (1988) observed large red asymmetry of the Ha line during 
the period of hard X-ray burst (Fig. [8J). These asymmetries were used to deter- 
mine the downward velocities estimated from the maximum shift of the centroid 
of the bisectors. The averaged over all red-shifted pixels during the impulsive 
phase velocity was of the order of 60 ± 10 kms . The downflow analysed in 
the Ha data and the upflow observed in the X-ray lines allows the authors to 
analyze the momentum balance of the flow. They conclude th at the momenta of 
upflow ing and downflowing plasma are approximately equal. IZarro &: Canfield 
(1989) conclude that the downflow velocity measured from red wing enhance- 
ment can be used as a diagnostics of impulsive solar flare heating conditions 
(Fig. ED. 

The work of iDing et al.l ( 19951 ) shows that the velocity of chromospheric 
downflows deduced from the red asymmetry of Ha line is around 30 - 40 kms -1 
with the lifetime of the order of 2 - 3 minutes. There are two major problems 
found by the authors: Why is the line center nearly not shifted while the line 
wing shows great asymmetries? The second problem concerns the life time of 
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Figure 8. The hard X-ray emission (left) and the Ha excess line profiles 
(right) observed during the impulsive phase of a solar flare. Dashed lines mark 
the centroids of the profile. The obtained downward velocities are around 
50kms _1 . The numbers on the X-axis of left panel represent minutes after 
23 U T and the arrows - the times of Ha observed line profiles (|Zarro et al.l 
Il988h . 




Figure 9. An example of the Ha line profiles observed during the impulsive 
phase of a solar flare. The red wing enhancement is observed at the same 
time as strong hard X-ray emission. Note th at the central part of the line is 
slightly blue-shifted (jZarro fc Canfieldl[T989l) . 



the do wnflow which is considerably longer than the life time predicted in lFisher 
(1989) simulations. Recent hydrodynamic and radiative transfer simulations 
may now explain these two problems. 

There are more papers presenting the observations of the red asymmetry 
related to the chro mospheric condensations driven during the impulsive phase 
of solar flares (e.g., IWuelser Sz MartJ ( 19891 ) - Fig. [TO]) . All of them determine 



more or less consistent observational picture of the chromospheric flows during 
the impulsive phase of flares. This picture is based on many spectroscopic ob- 
servations of the chromospheric line profiles. An important step was done when 
it became possible to calculate theoretical line profiles and compare them with 

obse rvations. 

ICanfield Gavlevl ( 19871 ) computed time-depen dent Ha line profiles for the 
dynamic model atmosphere of iFisher et al.l ( 1985al ). They simulate the effects 
of power-law electron beam heated chromosphere. Solving the radiative transfer 
equations for one-dimensional model atmosphere the evolution of Ha line profile 
was estimated. The time of the electron beam heating was 5 s and for detailed de- 
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Figure 10. The hard X-ray emission (top) and the Ha line profiles with the 
red asymmetry (bottom) observed during the imp ulsive phase of a solar fl are. 
Dotted lines mark the reference quiet Sun profile (Wu elser fc Martilll989f ). 



script ion of other parameters and computational methods see lCanfield &: Gavlev 
(jl987l ). In Fig. [□] the time sequence of the Ha line profile is presented. Dur- 
ing the non-thermal heating, the red-shifted component is present but after the 
heating was stopped, the Ha line exhibit the blue asymmetry although the cen- 
tral absorption feature is shifted towards longer wavelengths. This behaviour is 
explained by downflow of the chromospheric condensation. It is also worth to 
notice that the response of the Ha emission to the non-thermal electron beam 
is very fast (less than second). 

Similar but more precis e simulations of the dynam ics and radiation in a solar 
flare loop was presented by lAbbett Sz Hawlev ( 19991 ). Except the non-thermal 
heating of the chromosphere, they took into account the thermal hea ting by the 
soft X-ray irradiation within 1 - 250 A range. Carlsson Stein] ( 19971 ) radiative- 
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Figure 11. Time evolution of the calculated Ha line pr o files for the electron 
beam heated model atmosphere of iFisher et alJ (|1985a[ ) (jCanfield k, Gayleyl 
[1983). 



hydrodynamic code was used to analyze the response of the lower atmosphere 
at the footpoint of a flare loop. In the radiative transfer calculations the impor- 
tant transitions of hydrogen, helium and singly ionized calcium and magnesium 
were treated in non-LTE. One-dimensional atmospheric model was used in the 
calculations. 

As a starting models the authors took two different cases PF1 and PF2. The 
temperature and electron density stratifications of both preflare atmospheres are 
shown in Fig. [T2J Three levels of the non-thermal heating was considered which 
correspond to weak (F9), moderate (F10), and strong (Fll) non-thermal flare 
heating. The PF1 atmosphere is heated for 70 s with the F9 and F10 fluxes, 
and the PF2 atmosphere is heated for a shorter, 4 s burst but with strong Fll 
heating. 

Figure [T3l presents the time evolution of emergent Ha and CallK line pro- 
files. For the Ha line separated blue-shifted component is clearly visible while 
Call line exhibit red-shifted component. The contribution function calculated 
for these two lines explain why we observe such a two-components and asym- 
metric profiles (Fig. [Tl"|) . 

This analysis shows that the evolution of non-t hermally heated chromo - 
sphere progresses through two distinct dynamic phases (jAbbett fc Hawlevlll999l ) : 
a gentle phase, where the non-thermal energy input of the flare is essentially ra- 
diated away into space, and an explosive phase, where the flare energy rapidly 
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Figure 12. The temperature and electron density stratifications of both 
preflare atmospheres PF1 and PF2 compared to the standard semiempirical 
VAL3C chromospheric model of Vernazz a et al.l (119811) and to the semiemp ir- 
ical active atmosphere MF1ME of lMetcalJ (|199oh (|Abbett fc Hawlevi [T999) . 
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Figure 13. Left: Time evolution of the Ha and CallK line profiles for 
moderate heating of the atmosphere by non-thermal electrons (model F10). 
In each panel, the vertical axis represents relative intensity with respect to 
the continuum level and the horizontal axis denotes the wavelength from the 
line centre. The dashed lines in each panel represent the preflare line profile. 
Right: Time evolution of the Ha and Ca II K line profiles bu t calculated for 
strongly heated atmosphere (model Fll) ([Abbett fc Hawlevlll999f ). 



heats the atmosphere and drive large amounts of chromospheric material up 
into the corona, and down toward the photosphere. During the explosive phase, 
there is significant plasma motion and there are steep velocity gradients. More- 
over, the effects of thermal X-ray heating of the chromosphere remain negligible 

compared to the non-thermal heating in the impulsive phase. 

Similar, but much more extended calculations were presented bv lAllred et al.l 

( 20051). The basic s of co mputational methods are similar to those described by 



Abbett fc Hawlevl (jl999h but there are some significant improvements. The au- 



thors include the double power-law electron beam energy distributions recently 
observed in solar flares with the Reuven Ramaty High-Energy Solar Spectro- 
scopic Imager (RHESSI) satellite. Additionally, the effects of XEUV heating 
from a large number of high-temperature lines was taken into account using 
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Figure 14. The approximate formation height of different parts of the 
spectral line profile may be described by contribution function. This fig- 
ure presents components of the intensity contribution function for CallK 
(left) and Ha (right) lines after 50 s of flare heating in model F10. Line 
frequencies are in velocity units. Atmospheric velocities (dashed lines) are 
taken to be positive toward the corona. Thus, negative velocities associ- 
ated with downward-mo ving material correspond to red-shifts in the profile 
IjAbbett fc Hawlevl[l999h . 
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Figure 15. Left: Comparison of the electron beam heating rate in the pre- 
flare atmosphere for the F10 flare mode l. The solid line shows the heating rate 
used in the paper of lAllred et alJ (120051). and th e dashed line shows the initial 
beam heating rate of lAbbett fc Hawlevj[1999l ). Right: Solid line - thermal 
XEUV heating used in lAllred et all (l2005lT dash ed line - soft X-ray heating- 
rate used previously by lAbbett fc Hawlevl (|1999( ). 



results from the CHIANTI and ATOMDB databases and a wide range 1 - 
2500 A was used for direct thermal heating of the chromosphere. Figure [T5l 
pres ent the comp arison of the heating rates used by lAbbett &: Hawlevl ( 19991 ) 
and lAllred et all (|2005l ) . 

These new calculations confir med the previous r esults . However, the line 
profiles evolution differs from the lAbbett h Hawlevl ( 19991 ) calculations - the 
line asymmetry is not so significant and the blue and red components of the 
Ha and CallK lines, respectively, are not observed separately. Instead, the 
lines are asymme t ric w ith blue or red wing enhancement (Fig. [TBI) . As in 
Abbett Sz Hawlevl ( 19991 ). the authors found that the impulsive flare naturally 



divides into two phases, an initial gentle phase followed by a period of explosive 
increases in temperature, pressure and velocity. 
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Figure 16. Time evolution of the synthetic line profiles of Ly-a, Ha, 
Hell 304 A, and CallK lines for the moderate level of non-thermal heat- 
ing (model F10). Times are indicated at the top of each column. The dotted 
lines indicate the level of the continuu m close to the li ne center, while the 
dashed line is the preflare line profile (|Allred et al.ll2005l ). 



There are more papers which tr eat the problem of plasma flows i n the 



non-thermally heated ch r omosphere (Mariska et al. 1989; lEmslie et al.l Il998l : 



iKarlickv k Henouxl 1 19921 : iGan k Fand Il990l : iGan et al.1 Il99ll 1. All these hy 



drodynamic simulations predict an upflow of the hot coronal plasma due to the 
enhanced pressure in the region heated by non-thermal electrons or protons. 
This upflow is associated with downflow of chromospheric condensations, but 
with much lower velocities (Fig. I17p . These condensations disturb the line pro- 
files emitted from the chromosphere and cause significant asymmetries observed 
e.g., in Ha line. 

The calculated line profiles emerging from flaring atmosphere show roughly 
similar behaviour than the observed ones but unfortunately, the appearance 
and the time evolution of the calculated line profiles was not compared to the 
spectroscopic observations of particular flares. The validity of the modeling 
summarized in this chapter was not confirmed observationally up to now. 
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Figure 17. Left panels: Hydrodynamic response of the temperature and 
velocity to the non-thermal electron beam heating (iMariska et al.lll989f) an d 
to the non-thermal proton beam heating (right panels) (jEmslie et al . 1998). 
Positive values of the velocity correspond to upward motion of the plasma. 
Labels 10, 20, and 30 are the times in seconds after start of non-thermal 
heating. 



4. Velocity Field in Semiempirical Models of the Flare Atmosphere 

Asymmetries observed in chromospheric line profiles are also modeled using 
semiempirical flare models. This approach is based on the idea that the temper- 
ature stratification of the atmosphere is determined empirically in the way to 
reproduce the calculated spectrum in the best agreement with observations. This 
means that the energy-balance equation is not considered. Many spectral lines 
and continua are used to construct such kind of models. I nitially, static semiem - 
pirical models were developed for a quiet-Su n atmosphere (Vernazza et al.lll98ll ) 



and then for the static flaring atmosphere ( Machado et al. 1980T ). These mod- 
els were constructed under assumption of 1-dimensional geometry and hydro- 
static equilibrium (Fig. I18p . They are static however the time sequence of many 
semiempirical models can be used to describe the evolving atmosphere but this 
method is valid only for slowly-evolving atmospheres. More detailed description 
of semiempirical models can be found in the paper of Mauas in this book. 

Non-LTE radiative transfer methods applied to semiempirical models al- 
lows us to calculate the spectrum emerging from the atmosphere. In particular, 
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Figure 18. Temperature as a function of colum n mass for the flare mo dels 
Fl and F2, for the quiet-S un QS (VAL-C) mo del of lVernazza et alj (1 198 If) and 
for the plage (P) model of iBasri et alj (| 19791 ) (from iMachado et alj 1980). 



the p rofiles of the chrom e-spheric optically thick lines may be calculated in de- 
tails ( Heinzel et al.l fl~994j ) . It is also possible to reproduce asymmetric line pro- 
files solving the transfer equation for a dynamic atmosphere with velocity field. 



previously calculated level populations for static model (|Berlicki et al 



using 

2005). However, this approach is justified only for relatively small velocities 
(V < 10 km s -1 ) which do not significantly affect the level populations of the 
static model ( Nej ezchlebal 1 1 9981 ) . It cannot be used for impulsive phase of flares 
to model the chromospheric condensations which move quite fast. Therefore, 
such simplified calculations are used to model e.g., the gradual phase of solar 
flares, when the velocities in the chromosphere are low. 

One of the first semiempirical modeling of chromospheric flows was per- 
formed bv lGan et al.l ( 19931 ). Using the Ha line profiles observed for two flares 
the authors constructed the series of semiempirical models with chromospheric 
condensations. It was shown that chromospheric condensations are responsible 
not only for the red asymmetry of the Ha line, but also for the blue asymme- 
try of the line with central reversal (Fig. [T9j) . Chromospheric condensations 
were assumed to be homogeneous with constant velocity. The most important 
conclusion of this paper is that the properties of chromospheric condensations 
seem to be consistent with the results of hydro dynamical models of solar flares. 
Comparison of calculated Ha line profiles with real observations present also a 
valuable part of this paper. 

An interesting work was presented by iNeiezchlebal ( 19981 ) who simulate the 
influence of the velocity field on the Ha line profiles. The calculations were 
performed using non-LTE model of plane-parallel solar flare atmosphere with 
stationary velocity field. This velocity field was applied to different layers of the 
solar atmosphere and the e mergent Ha line profil es were calculated for two mod- 
els of solar flare Fl and F2 (|Machado et alJll980h . Figure 1201 present an example 



of asymmetric line profiles calculated for a weak-flare model atmosphere Fl and 
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Figure 19. Two examples of semiempirical models of the chromosphere with 
cool condensation. Left three panels correspond to the condensation which 
is responsible for the red asymmetry of the Ha line while right three panels 
presents condensation producing Ha line with blue asymmetry. Temperature, 
electron density and velocity stratification is presented for both cases (negative 
sign of the velocity corresponds to downflow). The calculate d Ha line profile s 
are plotted with solid line, observed ones - with dotted line ()Gan et al.lll993h . 
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Figure 20. Left: Models of the velocity field in the chromosphere used in 
the non-LTE calculations (positive sign of the velocity corresponds to down- 
flow). Right: Some examples of synthetic Ha line profiles calculated for a 
static and dynamic atmosphere with different velocities 10, 30, and 50kms _1 
(|Nejezchleballl998f ). 



for different velocity fields. The main conclusion of this work is that the velocity 
field affects the level populations via the increase of the downward radiation. 
Nevertheless, for velocities that do not exceed the thermal velocity of plasma, 
one can use the static populations for the formal solution of radiative trans- 
fer equation including the velocity to reproduce the observed line asymmetries. 
Other important point is that application of the bisector method would lead in 
some cases to reverse velocity, in others to underestimation of the velocity. It 
includes, besides the part of the profile directly affected by the moving material, 
also a "static" part of the profile. To use the bisector in terms of Doppler-shift 
the static part should be somehow eliminated. This remark makes question- 
able all estimations of the Doppler velocity obtained with the bisector method 
applied to self-reversed or emission chromospheric lines observed in solar flares. 
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13:44:00 13:44:30 



Figure 21. An example of the observed (thin line) and fitted (thick line) 
H(5 line profiles and determined velocity field across the atmosphere. Negative 
value indicates an upward velocity. In the lower panel the time evolution of the 
velocity value at two different heights in the atmosphere is presented. Filled 
circles refer to the hei ght « 900 and open circles to the height « 1400 km 
(jFalchi & Mauasll2T)0l . 



The conclusions of iNeiezchlebal (1998) suggest that deducing the velocity 
from flare line profiles is rather difficult and cannot be done only by searching 
for Doppler-shifts with the bisector method. Chromospheric line profiles suggest 
that the flare atmosphere is highly dynamic and stratified with rather compli- 
cated plasma motion. These lines are optically "thick" and the only reliable 
way to analyze the flows is to use the non-LTE radiative transfer codes, which 
enable us to compute the chromospheric models with velocity fields. Resulting 
synthetic line profiles can then be compared with the observed ones. 

The direct compa rison of the observed an d synthetic line profiles was pre- 
sented in the paper of Falchi Mauad (2002). They study the chromospheric 
structure of a small flare and construct 5 semiempirical models for different 
times, which reproduce the profiles of the H5, CallK, and Si I 3905 A lines dur- 
ing the flare evolution. In order to reproduce the asymmetry of the lines the 
velocity fields were introduced in the lin e profile calculations. T he modeling was 
done using the non- LTE Pandora code of lAvrett fc Loeserl (|1984| ). The trial-and- 
error method was used to reproduce the observed line profiles by the synthetic 
ones. Figure [2T1 (upper panels) present an example of the observed and fitted H5 
line profiles and determined velocity field across the atmosphere. In the lower 
panel of Fig. [21] the time evolution of the velocity value at two different heights 
in the atmosphere is presented. The presence of an upward motion in the flaring 
atmosphere at 1400 km, might be a signature of the chromospheric evaporation 
observed at chromospheric levels. It is interesting to notice that around 13:44:30 
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Figure 22. Left: Simulation of the influence of the downward velocity field 
on the emergent Ha line profile emitted from flaring model atmosphere Fl of 
Machado et al. (1980). In the middle panel the velocity used in the modeling 
as a function of temperature is shown (positive sign of the velocity indicates 
downflow) . Also an exa mple of the observed Ha line with red asymmetry is 
presented (for details see lHeinzel et al.l (|1994D V 



UT the downflow is observed at the height of 900 km, while the upward motion 
is evident at 1400 km above the photosphere. One possible explanation is that 
chromospheric evaporation together with condensations is observed. In this case 
chromospheric evaporation is observed at chromospheric levels and not, as more 

common, at coronal levels. 

Using th e semi empirical Fl model of a weak flare ( Machado et al.| [l980). 
Heinzel et al. ( 19941 ) showed that the blue asymmetry of Ha line profile is ob- 



served due to the downflow of chromospheric plasma. It is interesting to notice 
that the blue asymmetry is associated wit h the red-shift of the central absorp- 
tion feature. Similar results was shown by iGan et al. ( 1993 ) . The structu re of 
the velocity field used in the non- LTE simulations of lHeinzel et all (|l994l ) was 
qualitatively consistent with the concept of downward-moving chromospheric 
condensations (Fig. [22l). T h ese ca lculations were performed using the non- LTE 
code developed by iHeinzel (1995) and modified for flare conditions. The code 
uses a ID plane-parallel geometry and the atmosphere is in hydrostatic equilib- 
rium. Hydrogen excitation and ionization equilibrium have been computed by 
solving simultaneously the radiative transfer equations, the equations of statis- 
tical equilibrium for a five-level plus continuum atomic model of hydrogen and 
the equations of particle and charge conservation. The equations of statistica l 
equilibrium have been preconditioned according to iRvbicki Hummer! ( 199 ll ). 
The preconditioning is based on the lambda-operator splitting technique, where 
the exact lambda operator is expressed as an approximate operator plus the 
correction. Then the correction is iteratively applied to a lagged source func- 
tion by using the so-called Accelerated Lambda Iterations (ALI) method. For 
multilevel atoms thi s method is referred to a s MALI - Multilevel Accelerated 
Lambda Iterations ( Rvbicki Sz Hummer 199ll ). The preconditioned equations 
are then linearized with respe ct to the ato mic level populations and electron 

density and solved iteratively ( Heinzel 19951). 

This non- LTE code was also used in Berlicki et al. ( 20051 ) to analyse the time 
evolution of the line asymmetry observed during the gradual phase of the solar 
flare on October 22, 2002. In this paper for the first time the evaporative flows 
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Figure 23. The Ha line-centre optical depth distribution of the velocity field 
used in the modeling plotted for different values of r m and for Vq = —8 km s _1 
(upper left panel) and for different values of Vq and for r m = 1.0 (upper right 
panel). In the lower panels the influence of the velocity field on Ha line profiles 
emitted from flare is presented (upflow defined by Vo = — 8 and downflow 
Vq = +8 km s _1 for three values of l ogr m = 0.1, 1.0, and 10.0 (dotted, dashed 
and continuous lines, respectively) (Berlicki ct al. 2005). 



in the gradual phase are studied quantitatively by using a non-LTE radiative 
transfer code and spectroscopic observations of the flare ribbons. First the 
authors analyse the influence of different velocity fields on the emergent Ha line 
profile. Again, it was shown that the downflow of flaring plasma causes blue 
asymmetry of the self-reversed line while upflow - red-asymmetry (Fig. [23]). For 
the modeling of the observed line asymmetry except the changes of the value of 
the velocity, also the height of the velocity field in the atmosphere was different. 
The procedure of fitting the Ha line profiles was performed using a grid of many 
models by varying different parameters. Each observed profile was fitted by the 
least-square technique to a closest synthetic profile from the grid and the model 
with the velocity field was found for each analysed line profile. 

In the analysis the MSDP (Multichannel Subtractive Double Pass) spec- 
trograph (Mein 1991) coupled to the VTT telescope working at the Teide Ob- 
servatory (Tenerife, Canary Islands) was used. 36 Ha line profiles (six areas at 
six different times) observed during the M1.0 flare on October 22, 2002 were 
taken for the analysis (Fig. [25]- upper left panel). As an example we present in 
Fig. [25] (lower panel) some profiles observed in chosen area at three times. These 
observed profiles (solid lines) are fitted with the synthetic ones (dashed lines) 
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Fi gure 24. The imag e of the flare on October 22, 2002 used in the analysis 
of lBerlicki et a l. (2005) and the observed (continuous lines) and fitted (dashed 
lines) Ha line profiles (upper and lower left panels). 0-5 are the areas used in 
the analysis. The time evolution of the velocity in the chromosphere deduced 
from line asymmetries is plotted for two analysed areas in the right panel. 
Negative velocities correspond to upflow (|Berlicki et al.|[2005l) . 



obtained from the grid. In the right panel of Fig. [23] the temporal evolution of 
the velocity field is presented for two different areas of the flare. 

The authors inter p ret th e upflows found in the flare ribbons in terms of the 
Antiochos Sz Sturrockl ( 1978I ) model for gentle evaporation. This process may 



occur during the gradual phase of solar flares and it can be driven by conductive 
heat flux from the high-temperature flare plasma contained in magnetic flux 
tubes above the photosphere. In the future it would be interesting to use more 
spatial points at more times and to use the spectra obtained within a wider range 
of wavelengths. Other distributions of the velocity field in the chromosphere 
should also be tested. In addition, to perform non-LTE modeling of the flare 
structure it would be useful to have other spectral lines formed at different levels 
of the chromosphere. 



5. Summary 

In this review I presented some interesting papers concerning plasma flows ob- 
served during solar flares in cool chromospheric layers. These flows are directly 
responsible for the line-asymmetries and/or line-shifts often observed in chro- 
mospheric lines emitted by the flaring plasma. 

An important work was done to understand the flows and their mechanisms. 
In order to determine the plasma velocity and flow direction the bisector method 
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was applied for line profiles. Unfortunately, as we could see, this method leads 
in some cases to misleading estimations of the velocity. Recently, the direct 
comparison of the observed and synthetic line profiles gives more valuable in- 
formation about velocity fields in the chromosphere. All the data support the 
evaporative model of solar flares where explosive chromospheric evaporation of 
the hot plasma is associated with the chromospheric condensations observed in 
"cool" chromospheric lines. In the late phases of flares the gentle evaporation 
may be observed in chromospheric lines. 

For the future it is necessary to use large and dense grids of the chromo- 
spheric models computed with hydrodynamic and non-LTE codes. They may 
help us to understand the flows and give more realistic description of the phys- 
ical processes during the flares, particularly the heating mechanisms and their 
role at different phases of the flare evolution. 

Finally, really good spectral observations of flares are needed. They have 
to be co-spatial, simultaneous and obtained in different spectral ranges (X-ray, 
EUV, UV, optical, IR). Such observations would be very helpful to construct 
the full picture of the plasma flows during flares. There are some data concern- 
ing the flows observed in soft X-ray and EUV but they are extremely rare and 
almost never cospatial nor simultaneous with the observations in chromospheric 
lines. 
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